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1. Motivations 

Ultimately, we would like to study the finite temperature phase transition in QCD and the prop- 
erties of the plasma phase (QGP), with physical values of the quark masses and in the continuum 
limit. The present project attempts to set the stage for such a complete study. 

In this study we consider QCD with two degenerate flavors. A cross check of the lattice results 
obtained with staggered and Wilson fermion action will help controlling discretization artifacts, 
while a twisted mass term is expected to facilitate the continuum and chiral limits. 

Consider the phase diagram of two-flavor QCD in the temperature-mass plane: the first order 
deconfinement transition stemming from the infinite mass (or quenched) theory weakens with lower 
quark masses, until it turns into a smooth crossover for intermediate quark masses. In the chiral 
limit there has to be a true phase transition again, but its nature is still under investigation [[I]]. 

As a first step of our program, we wish to find the location of the phase boundary between 
hadronic and plasma phase, i.e. the pseudo-critical temperature and mass combinations, while 
taking advantage of the properties of twisted mass QCD. This means that we will have to explore a 
three-dimensional space of temperature T, bare quark mass m, and twisted mass parameter ju. 

2. Why Twisted Mass QCD Thermodynamics ? 

Wilson fermions have several advantages over staggered fermions, but they also have a more 
subtle chiral behavior, and a complicated phase structure, both at T = [|2[ and at finite temper- 
ature 

The twisted mass approach improves over the standard Wilson behavior in two 
ways: first, it prevents the occurrence of exceptional configurations and should make it relatively 
easy to reach mass values of the light pseudoscalar mesons close to the physical pion mass; second, 
once the Wilson hopping parameter fc is set to its critical value, the twisted mass term behaves as a 
conventional quark mass, and, at the same time, an 0(a) improvement is automatically guaranteed. 
For recent results see Refs. |M |9|] and for a review Ref. JT0|]. 

In this first report we search for the transitions between the hadronic and plasma regimes by 
varying the Wilson hopping parameter k related to the bare quark mass by K = 1/ (2m + 8) at fixed 
/3 and fixed twisted mass parameter }X. 

3. Strategy and simulations 

The simulations were performed on a 16 3 x 8 lattice with an improved version of the HMC 
algorithm as detailed in Ref. JTI| ] and with a Symanzik tree-level improved gauge action. They 
used approximatively three months x crate of apeNEXT jl2[]. We choose to work at j8 = 3.75 and 
/3 = 3.9 in order to take advantage of the T = results [||, ^]. In principle, we can then cross the 
(pseudo-)critical line at a fixed temperature by tuning the quark mass, either by varying K and/or 
IX. 

For this strategy to be successful, the simulation parameters N t and /3 need to satisfy: 

jchiral < simulation = { /(ty fl (j3)) < T ? Uenched . (3. 1) 

p or simulation > quenched ^ hadronic phase cannQt exist> while for simulation < ^chiral ±e QQp 

cannot exist, irrespective of the mass value. 
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Figure 1: (P) as a function of K (left diagram); Conjugate Gradient (CG) Iterations superimposed with 
(magnified) as a function of K at fixed /3 =3.75 and ji = 0.005 (right diagram). The results indicate a 
thermal transition or crossover at K t = 0. 165(1 ). 



We fixed N t = 8 by taking into account the lattice spacing from the T = studies, namely 
a(3.75) ~ 0.12 fm, a(3.9) ~ 0.095 fm, as well as the known critical temperatures, T° himl = 170 
MeVandr; rafW = 270MeV. 

In the first set of runs reported here we fixed ji. = 0.005 and varied the hopping parameter K. 



4. Results at j8 = 3.75 

At j3 = 3.75, the T = results show that the minimum pion mass which can be reached 
with our twisted mass parameter /J. = 0.005 is m n ~ 400MeV extrapolating existing results at 
/i = 0.005 [jSp. Our first goal here is merely to check whether a thermal phase transition or a 
crossover can be found in the required range K t < K c . 

Figure |l| (left) shows a scan of the average plaquette as a function of k. The steepest slope of 
the plaquette as well as the increase of the number of CG iterations needed for the inversion shown 
in Figure [I] (right), both suggest a crossover or phase transition at [JSj] 

jq(/3 = 3.75, jU = 0.005) = 0.165(1) . (4.1) 

Hence K t < K C (T = 0) = 0.1669, as required. 



5. Results at /3 = 3.9 

After the exploratory study at j8 = 3.75, the choice j8 = 3.9 brings us closer to the continuum 
limit. Results for T = at this j3 are available at a number of values for the twisted mass parameter 
jj., see Ref. [^]. The minimum pion mass at T = for our pL = 0.005, inferred from these results, 
is about 350 MeV. 

As a direct comparison between the results at the two /3 values we show in Figure 2 the number 
of Conjugate Gradient iterations required for the inversion - which is an indicator of criticality - 
as a function of K. 
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Figure 2: The number of Conjugate Gradient iterations as a function of K for the two j3 -values. The solid 
lines are a smooth interpolation to guide the eye. 



Figure 3 (left) shows a collection of results for the expectation value of the plaquette (P). 
The errors are smaller than the symbol, the solid line is a Bezier interpolation to guide the eye. 
We performed local fits to a straight line (P) = a + bK within subsequent intervals of width Afc = 
0.002, and we show in the same diagram the tangent in the inflection point. The parameters b are 
used as estimators of the derivative of the Plaquette w.r.t. to fc and are shown in Figure 3 (right). 
These results indicate a phase transition or crossover around k = 0.1597 located according to the 
maximal slope b max . To make this prediction more quantitative, our statistics was enhanced to 
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Figure 3: 

diagram) 



P) as a function of K (left diagram); -4^ as function of K at fixed /3 = 3.9 and fi = 0.005 (right 



0(10000) HMC trajectories on a selected sample of points: K = 0.1586,0.1591,0.1593,0.1597 
in the candidate critical region at j8 = 3.9. The results (Figure 4) suggest a long autocorrelation 
time in the critical region. Given these autocorrelations, our results are very preliminary. Most 
probably the errors are underestimated, but still the plots may serve as indicators for the location 
of a crossover or transition. 

Figure 5 shows the results for the Polyakov loop: the steepening is clearly seen, mostly thanks 
to the latest, high statistics results. The Polyakov loop increases in the same K interval as the one 
observed for the plaquette. The Polyakov loop histograms of the HMC results after thermalization 
are narrow in the two pure phases, and broaden around at K = 0.1597, indicating an increase of the 
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Figure 4: HMC evolution and error analysis for the high statistics runs at /3 = 3.9 : binned averages (left) 
as a function of the HMC trajectories, and results and errors as a function of the discarded HMC trajectories 
(right). 



fluctuations and a critical behavior. We see no two-state signal (two peaks in the histograms), which 
seems to exclude a first order transition, but only a finite size analysis can assess with confidence 
the nature of the critical behavior. The steepest slope of the plaquette and of the Polyakov loop, as 
well as the broadening of the probability distributions suggest a crossover or phase transition at K t : 

Kt(p =3.9, ii = 0.005) =0.1597(5) (5.1) 

and, as it was also observed at j3 = 3.75, 

^ < k c (T = 0) =0.16085 . (5.2) 

Although we postpone a detailed analysis of the spectrum and related observables to our on- 



going high statistics study Q13Q, we have collected a few results for the pion propagator. 
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Figure 5: The real part of the Polyakov loop as a function of K (left diagram), and the Polyakov loop 
histograms (right diagram) at j3 = 3.9 and fi = 0.005. The data set is the same as in Figure 2, with the 
inclusion of some high statistics results. Both diagrams are consistent with a critical point or crossover at 
K, =0.1597(5). 
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The (zero momentum) pion propagator G(t) is measured at a selected sample of couplings, 
and is fitted to a standard hyperbolic cosine form 



in the time interval [2:6]. 



G(t) =A cosM Af t 



N t 



(5.3) 
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Figure 6: The pion propagator G(?) for selected ff values. The solid lines are the simple fits described in the 
text. The right figure shows a subset of the same data points as in the left figure on a different scale. 

We show the quality of the results, with the fits themselves superimposed, in Fig. 6 (left and 
right, note a different scale between the two). The resulting fit parameters, A(mplitude) and M(ass), 
are collected in Fig. 7 . They indicate that the effective pion mass decreases while approaching the 
thermal transition from below, while the amplitude of the propagator increases, following the trend 
of the average plaquette. 

6. Summary and Outlook 

We have studied QCD with two flavors of dynamical Wilson fermions including a twisted mass 
term on a 16 3 x 8 lattice at two values of the temperature: j3 = 3.75 corresponding to T ~ 205 MeV 
and j6 = 3.9 corresponding to T ~ 259 MeV. 
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Figure 7: Amplitude of the propagator and effective pion mass from the hyperbolic cosine fits as a function 
of K. 
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In either cases we have simulated O(10) values of bare quark masses, by varying the hopping 
parameter K at constant pL = 0.005. We have observed a behavior consistent with a crossover (and 
not excluding a real transition) at a critical value of fc, which we denoted as K t , which is less than the 
critical K at T = 0. This behavior is similar to that observed with ordinary Wilson fermions [Q, ||]. 

In our current study Jl3| ] on the apeNEXT machines at DESY and INFN we have to study next 
the jit dependence of our results. To this end, we are repeating a K scan at a larger /J. value. At the 
same time we want to take full advantage of the twisted mass approach by working at full twist 
with k = K t .(j8, T = 0). In this case the phase transition will be crossed by tuning the twisted mass 
jU. Model studies along the lines of Ref. [Q] will be most useful to guide our simulations. 
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